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Spin is a welcome complication in the study of partonic structure that has led to new insights, 
even if theoretically and experimentally not all dust has settled, in particular on quark flavor de¬ 
pendence and gluon spin. At the same time it opened new questions on angular momentum and 
effects of transverse structure. In this talk the focus is on the role of the transverse momenta of 
partons. Like for collinear parton distribution functions (PDFs), we are also in the case of trans¬ 
verse momentum dependent (TMD) PDFs, talking about forward matrix elements. TMD PDFs 
(or in short TMDs) extend collinear PDFs with only spin-spin correlations to PDFs that include 
spin-momentum correlations, including also time-reversal-odd (T-odd) correlations, relevant for 
the description of single spin asymmetries. In this way TMDs open up new ways of studying the 
spin structure. Their operator structure within QCD, however, is more complex leading to various 
ways of breaking of universality. 
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1. Introduction 


Parton densities are a natural way of dealing with the notion of cross sections as an incoher¬ 
ent sum of scattering off the partons, quarks and gluons, in a hadron. They are not needed in 
electroweak cross sections for leptons where the targets and produced particles are the degrees of 
freedom of the lagrangian describing its interactions. Then one just needs the wave functions in¬ 
cluding spinors or polarization vectors to account for external states, plane waves including their 
possible polarization. At higher orders the calculations involves of course the whole regularization 
and renormalization procedures, but these can be handled in a renormalizable theory. 

The wave functions and spinors or polarizations of partons actually show up in the single pai‘- 
ticle matrix elements of the corresponding fields. For a proper treatment of parton densities, one 
needs to consider the matrix elements of the fields befween hadronic sfafes, e.g. a single hadronic 
fargef in inclusive deep inelastic scaffering (DIS), in fhe cross secfion leading fo forward mafrix ele- 
menfs of combinations of field operators. For local operafor combinations fhis one can sfill employ 
fhe full field fheorefical machinery. Such mafrix elemenfs, however, are merely momenfs of parfon 
densities. Helped by the kinematics in a high-energy process one can make a twist expansion by 
identifying the leading local operator combinations, which then can through a Mellin transform be 
identified wifh collinear parfon disfribufion functions (collinear PDF’s). The renormalizafion dif¬ 
fers for each of fhe local (composife) operators involving for each of fhem anomalous dimensions, 
leading fo mulfiplicative renormalizafion facfors for fhe momenfs and a convolution of PDFs with 
splitting functions to account for scale dependence of the matrix elements and the parton densities. 

This all works amazingly well for QCD, where the collinearity shows up in data, that exhibit 
jet structure linked to ’target remnants’ and ’scattered pardons’. There is a natural order of scales. 
Starting with the high energy scale of the process, say y/s and a fraction of it linked to the still 
high energy scale of the partonic process, say Q to the transverse momenta of the order of a GeV 
or less. Intermediate are larger values of the transverse momenta but viewed as part of the partonic 
probabilities these fall like as{\pW)/\pW- 

Of course high energy kinematics is an essential ingredient in this. The hadronic mass becomes 
irrelevant and a given hadron is struck at an instant (light-front time). If Pi • P2 —)• oo one has 5 —)• 00. 
Any momentum p = xiP{+X2P2 + Pt with pt-P\ = Pt'Pi = 0 and hadronic size p^ ^ pj ^ 
naturally must have xiX2 ~ M'^/s splitting into three types, xi finite and X2 ~ M'^/s —)■ 0 (p is 
momentum of a parton in 1), xi ~ /s —)■ 0 and X2 finite (p is momentum of parton in 2), or 
x\ ~ X2 ^Mj^/s —)• 0 {soft). For a parton in hadron 1 , one can then has a Sudakov decomposition 

p = x\PiP Pt p)^{p^ - p\)n, 

where n = P2/P\-P2- The approximate light-like vector n satisfies Pi-n = 1 andxi = p-n.\f anofher 
hadronic vector K for which P\-K s is used fo fix n, fhe change in x\ is of order js. The 
infrinsic fransverse momenfum pj will depend on fhe choice of n. Considering a process like fhe 
Drell-Yan process, fwo parfons of differenf hadrons combine info a hard final sfafe momenfum wifh 
(e.g. a lepton pair), the measurable quantity 


q = 


q-Pi 

P1P2 


Pi + 


q-P\ 

PyPi 


Pi + qT, 
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includes the experimentally accessible non-collinearity qr- The rest identifies the partonic momen¬ 
tum fractions (up to mass corrections) with scaling variables, 


x\ 


q-Pi 

PyPi 


IPi-q 


and X2 


q-P\ 
Pi Pi 


2 Prq' 


The basic subprocess is parton(pi) -i- parton (p2) —^ T{^)- The actual sensitivity to the transverse 
momentum shows up in the transverse momentum of the lepton pair. If p\ = x\P\ + p\T and 
P 2 = X2P2 + P 2 T one has qj = Pit + P2T at leading order linking the partonic momenta to the 
observed non-collinearity qj. The cross section in terms of the transverse momenta is a convolution 
in momentum space or a product in impact parameter space, for which we hope to write at leading 
order (in an expansion in M^/s) 


ria oc y d^piTd^P2TS^{PlT+P2T-qT)---^{xi,PlT',---)^{x2,P2T',---)--- 

oc J d'^br exp{iqT-bT)...^{xi-,bT-,...)^{x2-,bT-, (1.1) 


Besides hard amplitudes and kinematic factors the dots (in particular those in the arguments of the 
correlators) also include regulators and corresponding scale dependence to handle large pr UV 
behavior as well as rapidity divergences 0 . Effects of intrinsic transverse momenta of partons 
are best visible in (partially) polarised processes. In that case one has polarization vectors S for 
hadrons (parametrizing the spin density matrix) or measurable polarization vectors depending on 
final state distributions of decay products, e.g. in p —)■ Tin or in A —)• nN. These spin vectors can 
fix directions including in particular the transverse ones. The spin vector can be parametrized as 
MS = SlP+MSt —M^Stn, which obeys = —S\ + S\ and PA = 0 . Therefore the polarization is 
not correlated with the collinear component of the momentum. Collinear quark densities, indeed, 
just come as spin-spin densities, unpolarized densities (/f(x) or q{x)) in an unpolarized proton, 
longitudinally polarized (chiral) densities (gi(x) or ^q{x)) in longitudinally polarized proton and 
transversely polarized densities {h\{x) or 8 q{x)) in a transversely polarized proton. Similarly gluon 
densities are just unpolarized gluon densities (/f (v) or g(x)) in an unpolarized nucleon or circularly 
polarized gluon densities (gf (x) or Ag(v)) in a longitudinally polarized nucleon. 


2. TMD correlators and distribution functions 

The quark and gluon TMD correlators in terms of matrix elements of quark fields [|, ^ in¬ 
cluding the Wilson lines U needed for color gauge invariance of the TMD case are given by 

d>!y](x,pr;«) = f (2.1) 

2xrPP'^^^^{x,pr,n) = I^-^0^ (P,S|P"^(0)f/[o_^]P"^(^)f/f^ o, |P,S)|^^ (2.2) 

(color summation or color tracing implicit), where the Sudakov decomposition for the momentum 
pP of the produced quark or gluon is used. The non-locality in the integration is limited to the 
lightfront, ^-n = 0 , indicated with LF. The gauge links Ujq are path ordered exponentials needed 
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to make the correlator gauge invariant 0 ,^. Depending on the process under consideration differ¬ 
ent gauge links will appear [|0]- For the quark correlator the gauge link bridges the non-locality, 
which in the case of TMDs involves also transverse separation. The simplest ones are the future- 
and past-pointing staple links (or just [±]) that just connect the points 0 and ^ via lightcone 

plus or minus infinity, ±°°] ^[Or ^r] ^[±0° these as our basic building blocks. 

For gluons TMDs the most general structure involves two gauge links (triplet representation), de¬ 
noted as [t/,f 7 '], connecting the positions 0 and ^ in different ways. The simplest combinations 


More complicated possibilities, e.g. 


+] 


Uh 


- l! or its conjugate 


allowed for [U,U'] are [+,+’*'], [+,-’*'] and [-,-h^ 

with additional (traced) Wilson loops of the form f/Pl = U^ 0] ~ '^[0 ^ 

are allowed as well. A list with all type of contributions can be found in Ref. |j^ \fU = U' one 
can also use a single gauge link in the octet representation. 

Since the above correlator cannot be calculated from first principles, an expansion in terms of 
TMD PDFs is used, which at the level of leading twist contributions is given by [ jl^ |TI| , |I^ , lA] 






- {x, PT)r5$T + hiP {x, p\) ^ ^ ■ 


i 2' 


^PtSt 


-St'" ix,pl) + St'' i^^Pr) 


{x,Pt) + 


( U V 

PtPt mv_P 
m2 


2M2 


hf^^^\x,pl) 


-PtIpv} 


2M2 


O ^Is 


.pAfi^y} 


4 M 


— h\f\x,pl). 


( 2 . 3 ) 


( 2 . 4 ) 


We have used that +M^Sin^. For function like and the shorthand 

notation 

= Si,gfi]{x,pl) - ^^j^gfT{x,pl) ( 2 . 5 ) 

is used. The gauge link dependence in this parametrization is contained in the TMDs. Note that for 
quarks and are T-odd, while for gluons /j^^, hfj., and h^j are T-odd. 

Even if any gauge link defines a gauge invarianf correlafor, fhe relevanf gauge links fo be used 
in a given process jusf follows from a correcf resummation of all diagrams including the exchange 
of any number of A” gluons between the hadronic parts and the hard part, i.e. gluons with their 
polarization along the hadronic momentum. They nicely sum to the path-ordered exponential. For 
quark distributions in semi-inclusive deep inelastic scattering they resum into a future-pointing 
gauge link, in the Drell-Yan process they resum into a past-pointing gauge link, which is directly 
linked to the color flow in these processes. 


3. Operator analysis 

In the situation of collinear PDFs (integrated over transverse momenta), the non-locality is 
restricted to the lightcone, ^-n = = 0 (EC) and the staple links reduce to straight-line Wil¬ 

son lines. The correlators then involve the non-local operator combinations \j 7 { 0 )U^Q^^-^\l/{^)\LC or 
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f \ix:, expanded in terms of leading twist operators ■i/7(0)D” .. .D'’t/r(0) and 

Tr[f”^D'’.. .D"F"''(0)D” .. .D"] operators. Although gauge links are part of the matrix elements, 
they do not cause any non-universality or process dependence. In order to find the local operators 
for TMDs we integrate over transverse momentum including explicit transverse momentum vec¬ 
tors. Each transverse momentum becomes a derivative giving transverse indices. The simplest 
of these transverse moments are 


j = (3.1) 

Id^Pr (x) + cgi$g,,(x), (3.2) 

Id^pr 4>[^l(x,p,) = (x) +cgi^“^“^^^(x) (3-3) 


and similarly results for l(x,p 7 -). These integrated results will require standard UV regu¬ 
larization and corresponding scale dependence, while one also may need to consider appropriate 
combinations, e.g. subtraction of traces, to get finite results. Furthermore it can often be more 


appropriate to work with Bessel moments |14]. The correlators appearing in Eq. 3.3 are of the 
form 




d^-Pd^^T 

{2ny 




{P,S\Wj{0)U[0,^]O{^)Vi{^)\P,S) 


LF 


(3.4) 


where the 0{^) operators are rank two combinations of idj{^) = iDj[^) —Aj{^) and G“(i^), 
defined in a color gauge invarianf way (fhus including GEs), 


A?(^) = Jh-E£(^P-T].E)uW^jf«“(T])uW^,, (3.5) 

G“(^)= 


wifh £(Q being fhe sign funcfion. Nofe fhaf G“(i§) = does nof depend on ^ P, implying in 

momenfum space p-n = p^ =0, hence fhe name gluonic pole mafrix elemenfs |[T^, [I^, F/l, 18, T^, 
^0| ]. In Eq. ^]^one encounfers symmefrized producfs of fhese operafors indicafed wifh subscripfs 
{5G}, efc. Moreover fhe color summafion often infroduces mulfiple possibilifies, e.g. for a gluonic 
pole mafrix elemenf in combination wifh fwo gluon fields fhere are fwo possibilities, Tr(F[G,F]) 
(c = 1) and Tr(F{G,F}) (c = 2) fhaf have fo be summed over. The final imporfanf ingredienf in 
Eq. ^ are fhe gluonic pole facfors, calculable factors depending on fhe number of gluonic poles in 
fhe operafor and fhe pafh of fhe gauge link U. Mosf well-known are fhe single gluonic pole factors 

cW=±i. 


The operators ^q{x) can be considered as zero momenfum limifs or infegrafions over of mulfi- 
parfon correlafors such as fhree-parfon quark-gluon-quark correlators 4>o(x,y) wifh an operator 
sfrucfure (l^(0)D“(T])t/r(i^)), non-local along fhe lighfcone, and similarly ^p{x,y) or gluon-gluon- 
gluon correlafors like Ef (x,y). These mulfi-parfon disfribufions play a role in higher fwisf confri- 
bufions to fhe cross sections and can be used to esfablish relations, fhey appear in sum rules, fhey 
exhibif symmefries befween correlators involving parfons and anfi-parfons, efc. Several confribu- 
fions on fhese fopics are presented af fhis conference. 
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In many cases the multi-parton distributions are also driving the large pr-behavior of the 
TMDs. For the collinear PDFs extended to the corresponding TMDs, this is the a ^/pj behavior of 
/i {x,Py), in essence the evolution equations. In this way multi-parton distributions are driving the 


large pj behavior of many of the TMDs ||^, in particular the T-odd ones. Note, however, that 
the large pj behavor of TMDs without a collinear counterpart may be driven by collinear functions. 


such as /f driving the linearly polarized gluon distribution [ |23| , |24| ] or the transversity h 
driving the Pretzelocity function h\j [^1|]. 


4. Universal TMDs 


Taking transverse derivatives gives the coefficients in the expansion in transverse momenta, 
for which we like to use traceless iiTeducible tensors of a fixed rank, which describe in 

essence fhe azimuthal dependence. We would like to use the moments to identify the coefficients 
in the azimuthal expansion 


^{x,pt) = ^{x,p\) + ^ + ■ 


M 


(4.1) 


If in the higher moments just operators of the type <j>“i would appear, it would be easy to find 
the operator expressions for fi{x,p\). It would correspond with the rank zero operator ^{x,p\) 
including all d-d traces that are subtracted in the higher rank operators and account for the p\ 
dependence. Such is actually the case for fragmentation functions where the gluonic pole matrix 
elements (after integration over transverse momenta) vanish. 

For the distribution correlators, however, this procedure does not lead to a unique correlator 
linked to a particular function because gluonic pole matrix elements do not vanish, in particular also 
terms including G-G traces. For instance even the unpolarized quark (or gluon) TMD distributions 
f\{x,p^) remain gaugelink-dependent [^] because of this. Only the d-d traces are taken care 
of in the -dependence of the function, but gluonic pole trace operators d>G G,c(-^)Pr) need to 
be included and require introduction of functions 8/f^'^\x,p\). These functions have to satisfy 
/ d^pT Px df\^^’‘^\x,Px) = 0, so they cause gauge-link dependent modulations, 

(x,p\)= f\{x,p\)+ £ {x,p\) + ..., (4.2) 

c=1.2 

of which only the first term survives in the collinear, pr-integrated, situation. It leads to process 
dependence in the p\ dependence, e.g. in the p^-width. In this case there are two possible color 
contractions in the summation over c. There is no term with a single gluonic pole matrix element 
because of the time reversal nature of the gluonic poles being odd. The time reversal nature, 
however, makes T-odd functions gaugelink-dependent already from the start, such as for the Sivers 
function, 

{x,p\) = cf (x,p2 ) + .... (4.3) 

For the Pretzelocity distribution there are three different operator structures contributing to the 
gauge link dependence. 


hiP {x,Pt)= (x, pI)+ Y, ^GG,c{x,pl) + . 


(4.4) 


c=l,2 
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Thus measurements of Pretzelocity effects are gaugelink-dependent, and hence process dependent, 
even if the observable is a T-even function. In any given process a particular combination of the 
universal functions on the righthandside appears. These universal functions are here labeled by 
a combination of d and G identifying the operator structure and if needed an index c if multiple 
color configurations have to be considered. In all of the above expressions, one has to be aware 
of additional modulations that come from operators with even more (traced) gluonic pole terms. 
To study their possible importance lattice studies using different gaugelink structures would be 
useful 


30|. 


Turning to the gluon distributions, one also encounters gaugelink-dependence already for the 
unpolarized TMD distributions. 


C=1 


(4.5) 


including four different color configurations in the pj modulation. Like for the Pretzelocity, one 
has for a T-even situation such as that of linearly polarized gluons in an unpolarized hadron also 
multiple universal functions. 




(4.6) 


The above examples illustrate our ongoing efforts [Q] to establish a universal set of TMD functions. 
As a final note, I wanf fo sfress fhaf in sifuafions where fwo (or more) TMDs wifh nonzero rank are 
involved, one musf accounf for possible addifional color factors in fhe basic expressions fhaf are in 
DY-like processes for insfance differenf from fhe 1/Ac or 1/ (A/ — 1) factors for qq or gg inifiafed 
processes [| 


271. 


5. Conclusions 

TMDs encode many feafures fhaf can be linked fo fhe parfonic sfrucfure of hadrons and fhaf 
can in principle be accessed af leading order provided fhaf one picks fhe righf variable, usually 
involving azimufhal asymmefries in polarized processes. Such efforfs are under invesfigafion in fhe 
experimenfal programs af RHIC/Brookhaven, JLab, BELLE, COMPASS/CERN, JPARC, BESIII 
or BaBar (see confribufions af fhis workshop). In parallel, fheorefical developmenfs are underway 
fo undersfand fhe dafa and fhe way fhese have fo be incorporafed info our view of fhe parfonic 
sfrucfure of hadrons. This is a nonfrivial enferprise since fhere are many ideas, buf also many 
technical hurdles fo fake. In my falk I have focussed for a large parf of efforfs in esfablishing a 
universal sef, connected fo specific operafors, for TMDs fhaf one can fhen can fry fo work wifh [^]. 
I have nof addressed fhe extensive efforfs fhaf are ongoing fo undersfand fhe scale dependence and 
fhe mafching fhaf is needed fo simulfaneously undersfand fhe behavior af low and high qj values. 
Eurfhermore, links exisf wifh work fo undersfand low-x behavior, use of more involved hadronic 
observables like di-hadron fragmenfafion funcfions or double-parfon disfribufions fo undersfand 
mulfi-parfon processes, where again we refer fo ofher confribufions af fhis conference. 


7 





Spin Physics and Transverse Structure 


RJ. MULDERS 


Acknowledgments 

This research is part of the research program of the “Stichting voor Fundamenteel Onderzoek 
der Materie (FOM)”, which is financially supported by the “Nederlandse Organisatie voor Weten- 
schappelijk Onderzoek (NWO)” and the EU "Ideas" programme QWORK (contract 320389). 

References 

[1] J. Collins, “Foundations of perturbative QCD,” Cambridge University Press (2011). 

[2] J. C. Collins and D. E. Soper, Nucl. Phys. B194 (1982) 445. 

[3] J. C. Collins and D. E. Soper, Nucl. Phys. B193 (1981) 381. [Erratum-ibid. B 213, 545 (1983)]. 

[4] A. V. Belitsky, X. Ji and F. Yuan, Nucl. Phys. B656 (2003) 165. 

[5] D. Boer, P. J. Mulders and F. Pijlman, Nucl. Phys. B667 (2003) 201. 

[6] C.J. Bomhof, PJ. Mulders and F. Pijlman, Eur.Phys.J. C47 (2006) 147. 

[7] C.J. Bomhof and PJ. Mulders, Nucl.Phys. B795 (2008) 409. 

[8] M.G.A. Buffing, A. Mukherjee and P.J. Mulders, Phys.Rev. D86 (2012) 074030. 

[9] M.G.A. Buffing, A. Mukherjee and P.J. Mulders, Phys. Rev. D88 (2013) 054027. 

[10] P. J. Mulders and R. D. Tangerman, Nucl. Phys. B461 (1996) 197 [Erratum-ibid. B484 (1997) 538]. 

[11] A. Bacchetta, M. Diehl, K. Goeke, A. Metz, P. J. Mulders and M. Schlegel, JHEP 0702 (2007) 093. 

[12] P.J. Mulders and J. Rodrigues, Phys.Rev. D63 (2001) 094021. 

[13] S. Meissner, A. Metz and K. Goeke, Phys.Rev. D76 (2007) 034002. 

[14] D. Boer, L. Gamberg, B. Musch and A. Prokudin, JHEP 1110 (2011) 021. 

[15] A.V. Efremov and O.V. Teryaev, Sov.J.Nucl.Phys. 36 (1982) 140. 

[16] A.V. Efremov and O.V. Teryaev, Phys.Lett. B150 (1985) 383. 

[17] J-W. Qiu and G.F. Sterman, Phys.Rev.Lett. 67 (1991) 2264. 

[18] J-W. Qiu and G.F. Sterman, Nucl.Phys. B378 (1992) 52. 

[19] J-W. Qiu and G.F. Sterman, Phys.Rev. D59 (1998) 014004. 

[20] Y. Kanazawa and Y. Koike, Phys.Lett. B478 (2000) 121. 

[21] A. Bacchetta, D. Boer, M. Diehl and P.J. Mulders, JHEP 08 (2008) 023. 

[22] M. G. Echevarria, T. Kasemets, P. J. Mulders and C. Pisano, arXiv: 1502.05354 [hep-ph]. 

[23] S. Catani and M. Grazzini, Nucl. Phys. B 845 (2011) 291. 

[24] S. Catani and M. Grazzini, Eur. Phys. J. C 72 (2012) 2013, [Erratum-ibid. C 72 (2012) 2132]. 

[25] D. Boer, M.G.A. Buffing and P.J. Mulders, arXiv: 1503.03760 [hep-ph]. 

[26] M.G.A. Buffing and P.J. Mulders, JHEP 1107 (2011) 065. 

[27] M.G.A. Buffing and P.J. Mulders, Phys. Rev. Lett. 112 (2014) 092002 

[28] F. Hautmann et. al., Eur. Phys. J. C74 (2014) 3220, 

[29] B. U. Musch, P. Hagler, M. Engelhardt, J. W. Negele and A. Schafer, Phys. Rev. D85 (2012) 094510. 

[30] M. Engelhardt, B. Musch, P. Hagler, J. Negele and A. Schafer, PoS LATTICE 2013 (2014) 284. 



